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Abstract 

We investigate the Liouvillian integrability of Hamiltonian systems describing a universe 
filled with a scalar field (possibly complex). The tool used is the differential Galois group 
approach, as introduced by Morales-Ruiz and Ramis. The main result is that the generic 
systems with minimal coupling are non-integrable, although there still exist some values of 
parameters for which integrability remains undecided; the conformally coupled systems are 
only integrable in four known cases. We also draw a connection with the chaos present in 
such cosmological models, and the issues of the integrability restricted to the real domain. 

1 Introduction 

Homogeneous and isotropic cosmological models, although very simple, explain the recent 
observational data very well \57\ I55| . Their foundation is the Friedmann-Robertson- Walker 
(FRW) universe, described by the metric 
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where a is the scale factor, 6.^0.2 is the line element on a two-sphere, and we chose to use the 



conformal time rj. As can be seen from the above metric, the scale factor represents the relative 
change in the distance of two points whose spatial coordinates are fixed. It depends only on time, 
so that the whole universe is deformed in a homogeneous fashion. From the anthropocentric 
point of view it could be seen as a three dimensional space evolving in the time - in the simplest 
case when the curvature index K is zero, it would be a Euclidean space stretched according to a. 

If we were to fill such a universe with a matter, its properties could only depend on time, 
and be the same in all points of the spatial subspace at a given value of r/ - otherwise the model 
would no longer be homogeneous. 

For example a perfect fluid would be completely described by two quantities - its density 
and pressure as functions of time. A scalar fleld would be described by its field variables. Also 
a cosmological constant with its trivial dependence on time could always be included in such 
models. 

Depending on the matter components one obtains various evolutions of the scale factor a, 
as given by the general action 



where TZ is the Ricci scalar, A the cosmological constant, V the field's potential, ^ the coupling 
constant, and g is the density of the perfect fluid. The potential usually includes at least a 
quadratic term m^|'0p, where m is the so-called mass of the field. When ^ = we say that the 
field is minimally coupled - it does not uncouple since the determinant of the metric g multiplies 
the whole Lagrangian density. Case with ^ = | is the so-called conformal coupling. 

For the considered geometry, the above action can be simplified so that it allows the Hamilto- 
nian approach with the phase variables depending only on conformal time r]. Due to the required 
covariance of general relativity, the system is subject to constraints, which in our case mean that 
the obtained Hamiltonian's value is zero. However, we note that including an additional matter 
component g is equivalent to considering other energy levels. Namely for g oc (which is the 
case for radiation) a constant is added to the Hamiltonian, thus imitating its non-zero value. 
This is the justification for studying the systems integrability on a generic energy hyper-surface. 

From the observational point of view, the cosmological constant provides an explanation 
for the current accelerating expansion of the universe (A Cold Dark Matter model [21]), but a 
better solution still is sought for. A real scalar field dubbed "quintessence" with the so-called 
slow rolling potential, which models the dark energy component has been extensively used for 
that purpose [lll[68|. Realisations of the field itself include also Bose-Einstein condensate of 
axions [25] or a phantom violating the energy principle '15]. 

Finally, scalar field could also be the mechanism behind the infiation |43^ I42j . which is 
currently the most established and used scenario for the early Universe [13] . Recently, Komatsu 
et al. [39] have shown that latest observational data (WMAP, SNIa, Barion Oscillation Peak and 
others) show that the model of chaotic infiation (which strictly speaking should be called non- 
integrable or complex in the sense that we demonstrate in the present work) with the quadratic 
potential remains a good fit (within the 95% confidence domain). 

From the physical point of view, a universe filled with only one component seems simple 
enough but it is not the case here. Chaotic scattering has been found in minimally coupled 
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fields [23j, as has chaotic dynamics |48j . 

The first of our results is that minimally coupled fields are not integrable in the generic case. 
There are however special families of the system's parameters which leave the question open. 
We give the appropriate conditions in the concluding section. 

There are several physical reasons to study more than just minimally coupled fields. Early 
works on chaotic inflation found the coupling constant ^ small or negative [30j but some argue 
[26) that the paradigm of inflation should be generalised to the case with non-zero coupling 
constant which should not be fine-tuned close to zero, and WMAP observations seem to indicate 
non- negligible 

The coupling could be generated by quantum corrections [101 [29j . or from the renormalisation 
of the Klein-Gordon equation as described in [16|. The coupling constant should be fixed by 
particle physics of the matter composing the scalar field, for example the way = 1/6 was found 
in the large approximation to the Nambu- Jona-Lasimo model in [3l] . Non-minimally coupled 
fields are also interesting in the context of description of the dark energy for which the ratio 
between the pressure and the energy density is less than —1. Such matter is called a phantom 
matter, and cannot be achieved by standard scalar fields [27] . 

Conformally coupled fields were subject to more rigorous integrability analysis, as opposed 
to minimally coupled ones, thanks to the natural form of their Hamiltonian. As will be shown 
in the next section, the kinetic part is of natural form, albeit indefinite, and the potential is 
polynomial (in the case of real fields). 

Chaos has been studied in such fields by means of Lyapunov exponents, perturbative ap- 
proach, breaking up of the KAM tori \17\ lllj. Also the Painleve property |i35] was employed as 
an indicator of the system's integrability. 

Ziglin proved that the system given by (j20p is not meromorphically integrable when A = A = 
and k = 1 [67]. His methods were also used by Yoshida to homogeneous potentials which is 
the case for the system when k = [611 [621 [631 [M] • Later, Yoshida's results were sharpened by 
Morales- Ruiz and Ramis |49] , and used by the present authors in [H] to obtain countable families 
of possibly integrable cases with some restrictions on A and A. Also recently, more conditions 
for integrability have been given in [T2], although only for a non- zero spatial curvature k and a 
generic value of energy, that is, when the particular solution is a non-degenerate elliptic function 
defined on a non-zero energy level. 

Our work shows, that the conjecture of that paper is in fact correct - as shown in Section 
4 - the conformal system is only integrable in two cases (with the above assumptions). We go 
further than that and show that for a generic energy value, a spatially flat {k = 0) the universe 
is only integrable in four cases. Also, the particular case of zero energy is analysed and new, 
simple conditions on the model parameters are found. Finally, we check that when E = k = 0, 
the problem remains open, as the necessary conditions for integrability are fulfilled. 

When it comes to numerical studies of the problem, there are various results, most notably 
a chaotic behaviour [36], but also a fractal structure and chaotic scattering [59]. However, it 
remains unclear whether the widely exercised complex rotation of the variables changes the 
system's integrability. Even for very simple systems it was shown ^31j that there might exist 
smooth integrals, which are not real-analytic. This question is especially vital since our Universe 
clearly has real initial conditions and dynamics. 

The results obtained for the conformal coupling are much stronger than in the minimal one. 
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We manage to show, that the four cases with known first integrals are the only integrable ones 
for the generic energy value. 

The Hamiltonian of both these systems has indefinite kinetic energy part, and to cast it 
into a positive-definite form a transition into imaginary variables is used. It has been done for 
a conformally coupled field [36], but some authors, see e.g. [51], argue that there are physical 
limitations which forbid extending the solutions through singularities such as a = 0, and an 
imaginary scale factor seems even less realistic. 

We would like to stress that the complexification of the variables in our approach is not to be 
connected with the physical evolution of the system into the imaginary values. The behaviour 
of differential equations in the complex domain is a tool that allows for obtaining general results 
regarding its integrability (also real), as can be seen on the example of the Painleve analysis 

Despite the fact, that systems of the considered type were often called non- integrable, there 
was seldom a rigorous proof of that proposition. However, the Liouvillian integrability can be 
studied successfully, as we try to show in this article. 

The authors are aware of one notable attempt at studying the problem in [21]; sadly, that 
paper contains a serious mistake in the application of Yoshida's theorem. The method used 
requires rescaling by the energy of the system, which the authors of [21] take to be zero. Thus, 
their Theorem 3 (which contradicts one of the results presented here) is in fact false. The zero 
energy level usually requires a separate treatment, which we also present here. 

The plan of the paper is the following. In the next section we describe the mentioned two 
cosmological models: of minimally coupled field and conformally coupled field starting from a 
general description up to the Hamiltonian formulation of these models. Section [3] is devoted 
to introduction to the Morales-Ramis theory. This is our tool to prove the announced rigorous 
integrability results for these models. In the next two Sections H] and [5] all details of integrability 
analysis are presented. For convenience of readers all integrability results are recapitulated in 
Section [6l 



2 Physical system's setup 

2.1 Minimally coupled field 

The general action ([2]) now includes the following parts 

1 



I 



WirG 



7^-2A 



m 



Using the coordinates of ([T|), the Lagrangian becomes 



£ = 6(aa + i^a2) + i|^|V-^a^|V'P 



-gd X, 



2Aa^ 



(3) 



(4) 



with the dot denoting the derivative with respect to time. We also dropped a coefficient which 
includes some physical constants and the part of the action related to the spatial integration. 
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Next we subtract a full derivative 6(da)', and use the polar parametrisation for the scalar 
field ip = -v/T2 </>exp(i^) to get 

2 

C = 6{Ka'^ - a2) + 60^(02 + ^2^2) _ ^IILa^^p-^ _ 2Aa\ (5) 
and obtain the Hamiltonian 

with 

Pa = -12d, p^ = 12a20, pe = 12a^4>'^e. (7) 

Note that the the elliptic constraints of general relativity require that the above Hamiltonian 
is identically zero - this is the so-called Friedmann equation, although it is not a dynamical 
evolution equation but rather a conservation law [47j . 

Since is a cyclical variable, the corresponding momentum is conserved so we substitute 
Pq = 2u;2. To make all the quantities dimensionless, we make the following rescaling 

m^^m^h^\K\, A^3A\K\, uj^ ^ 72u;^\K\, pl^72pl\K\, pi ^ 72pl\K\. (8) 

There is no need of changing the variables a and (j) along with their momenta, as this is really 
changing the time variable r], and thus the derivatives to which the momenta are proportional. 
This results also in dividing the whole Hamiltonian by 6\/2|-f^| to yield 

V2H = \ {-pI + Ip^) - A«2 + ^^4 ^ ^2^2„4 ^ 

If the spatial curvature is zero, any of the other dimensional constants can be used for this 
purpose, so without the loss of generality we take the right-hand side to be the new Hamiltonian 

H = \ I^-pI + ^pI^ - ka^ + Aa^ + m^'a^ + (10) 

and in all physical cases, k G {—1,0, 1}, ^ 0, G M, A G M, = 0. We extend the analysis 
somewhat assuming that the Hamiltonian might be equal to some non-zero constant E £M. We 
will later see, that our analysis includes also the possibility of these coefficients being complex. 

Note that for a massless (m = 0) field, the system is already solvable, as shown in Ap- 
pendix A. 

From this point on, we take u = 0, which means the phase is constant. Since the model has 
U(l) symmetry, we can always make such a field real with a rotation in the complex ■0 plane. 
In other words, we will be investigating a real scalar field only. The reason why we restrict the 
problem is the following: the method employed requires an explicit (non-constant) particular 
solution, and the only one known requires (p = 0; which is a singularity of the full Hamiltonian. 

Under the above assumption the Hamilton's equations of system (llOp are 

a = -Pa, Pa = 2ka - Aa'^ {A + 711^ <p'^) + \pI, 

.1 ^ (11) 

4> = —P<t>i P<P = -2m^a^<j). 
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We note that there is an obvious particular solution, which describes an empty universe: 
4' = P(p = = Pa = —Q- Thanks to the energy integral E = ^q^ + kq^ — Ag^, it can be 
identified with an appropriate elliptic function. 

2.2 Conformally coupled scalar fields 

The procedure of obtaining the Hamiltonian is the same as in the case of minimally coupled 
fields, only this time the action is 



167rG 



IV'r 



-gd X, 



(12) 



where an additional coupling to gravity through the Ricci scalar TZ, and a quartic potential term 
with constant A are present, as opposed to the minimal scenario. We keep the same notation as 
before and express the involved quantities in the same coordinates and get 

C = 6(aa + Ka') - \da\^\' + \\i;\'a' - ^a^P - ^a^l" - ^Aa^ - \Ka'\^\\ (13) 



from which we remove a full derivative, and introduce new field variables tp = \/i2 <?!'exp(z0)/a 
to obtain 



£ = 6 



2 « 

+ 0)0 —a +K{a — ) — tit a 6 — —a — Xo 

3 



The associated Hamiltonian is 



(14) 



H 



24 



P4> + laPe-Pa +6 



i^(</.2-a2) + ^aV + A0^ + |a^ 



(15) 



with 

p„ = -12d, p^ = 120, pg = l2<j?e. (16) 

We can see that is a cyclical variable because we took the potential to depend on the modulus 
of ip only, so we write a constant instead of the respective momentum pQ = uj. 

Finally, we express everything in dimensionless quantities, rescaling the constants, but also 
the time and momenta (as they are in fact time derivatives), which results in rescaling the whole 
Hamiltonian. We do this as follows 



A^^A|K|, \^]-X\K\, p1^1UpI\K\, H 



1 



12a/WI 



(17) 



when K ^ 0, and using another of the dimensional constants otherwise. Thus, eliminating a 
multiplicative constant, the Hamiltonian reads 



H = - 



I 



kicjP' - a^) + ^ + 171^0^4^ 



+ i (Aa^ + X<p^) 



(18) 



with k G {—1,0, 1} {K = k\K\); lo, X, A, G M, and H = in any physically possible setup. 
Exactly as in the previous case, the zero value of the energy is a consequence of the constraints 
introduced by general relativity. 
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We note that for m = the system decouples, and is triviahy integrable as shown in Ap- 
pendix B. That is why we will assume m 7^ henceforth. We will also take u; = 0, that is, 
consider a scalar field equivalent to a real field after a unitary rotation in the complex ip plane. 

We change the field variables into the standard q and p ones for further computation, taking 

a = qi, Pa=Pl, , . 



The Hamiltonian is then 



1 



h = -[-pI + pI)+v, 

\ 1 (20) 



3 Differential Galois obstructions to integrability 

Let (M, u) be a 2n-dimensional complex analytic symplectic manifold. For a meromorphic 
function H : M ^ C, we denote by Vh the Hamiltonian vector field generated by H and let us 
consider Hamiltonian equations 

dx 

— = vh{x), teC, xe M. (21) 

We assume that a non-constant particular solution ip{t) of system (I2ip is known. Its maximal 
analytic continuation defines a Riemann surface T with the local coordinate t. 

Linearisation of (I2ip around ip{t) yields variational equations of the following form 

^ = Am Ait) = ^iip{t)), eeTrM. (22) 

Thanks to Hamiltonian character of the system the dimension of variational equations can 
be reduced by two. First we use the fact that a Hamiltonian system has at least one first 
integral namely Hamiltonian H, thus we can restrict system (I2ip to the manifold = {x G 
M I H(x) = e}, where e = H{ip{t)). Then we consider the induced system on the normal bundle 
N := TvMe/TV of F 

7) = I(t)r/, A{t)ri = ^,{T{v){^-^i)), r] G N. (23) 

Here vr : T^M^ N is the projection. The system of 2n — 2 equations obtained in this way is 
called the normal variational equations. 

We can consider the entries of matrices A and A as elements of field /C := of mero- 

morphic functions on F. This field with differentiation with respect to t as a derivation is a 
differential field. Only constant functions from IC have a vanishing derivative, so the subfield of 
constants of /C is C. 

It is obvious that solutions of (I22p are not necessarily elements of /C". The fundamental 
theorem of the differential Galois theory guarantees that there exists a differential field J- IC 
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such that it contains n hnearly independent (over C) solutions of (I22p . The smallest differential 
extension ^ D /C with this property is called the Picard-Vessiot extension. A group Q of 
differential automorphisms of which does not change /C is called the differential Galois group 
of equation (j22p . It can be shown that ^ is a linear algebraic group. Thus, it is a union of disjoint 
connected components. One of them containing the identity is called the identity component of 

g. 

Differential Galois theory was created as a tool to answer the question: whether a given 
system of linear equations possesses a solution that can be written in a closed form, i.e. is it 
solvable? The main theorem of this theory states that the necessary condition of solvability in 
the class of Liouvillian functions (i.e. by generalised quadratures) is solvability of its differential 
Galois group. We can try to connect the integrability of the original nonlinear system with 
solvability of its variational equations. However, there is a more direct connection. Namely in 
eighties of XX century Ziglin observed that if system ()2ip has k > 2 functionally independent 
meromorphic first integrals, then variational equations ()22|) and also normal variational equations 
(|23p possess k rational first integrals and moreover the monodromy group (that is a subgroup of 
differential Galois group) has the same number of invariants [651 166j. Fourteen years later the 
relation between first integrals and invariants of the differential Galois group was analysed by 
Baider, Churchill, Rod and Singer in [19J. However the final formulation of relations between 
integrability of Hamiltonian systems and properties of the differential Galois group of variational 
equations due to Morales and Ramis |5H HS] where in their analysis not only the presence of 
first integrals is taken into account but also the consequences of the involution of first integrals. 
Their main theorem that will be the crucial tool of our analysis is the following. 

Theorem 1 (Morales- Ruiz and Ramis [49| ) Assume that a Hamiltonian system is mero- 
morphically integrable in the Liouville sense in a neighbourhood of a phase curve T and irregular 
singularities of the variational equations along F do not correspond to phase points at the in- 
finity. Then the identity component of the differential Galois group of the (normal) variational 
equations associated with T is Abelian. 

Let us explain assumptions concerning variational equations in the above theorem. Usually 
the Riemann surface corresponding to the phase curve T is not compact so we compactify it 
adding some points. Typically these points correspond to equilibria or infinite points. In the 
later case we have to add these points to the phase space, i.e, we have to extend our original 
system into a 'bigger' phase space. For the extended system the requirement that the considered 
first integrals are meromorphic in a neighbourhood of the phase curve put strong restrictions: 
they have to be meromorphic at the infinity. Thus if we remove the assumption about irregular 
singular points we have to restrict the class of first integrals. Below we give a version of the 
Morales-Ramis theorem without assumptions concerning the regularity of variational equations, 
which is adapted to Hamiltonian systems considered in this paper. 

Theorem 2 (Morales- Ruiz and Ramis [49]) Assume that a Hamiltonian system defined in 
a linear symplectic space is generated by a rational Hamiltonian function and is rationally in- 
tegrable in the Liouville sense. Then the identity component of the differential Galois group of 
the (normal) variational equations associated with F is Abelian. 
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In applications the most difficult part is to check the abelianity of variational equations. 
Fortunately, thanks to the separation of variational equations into two parts we can restrict to 
its normal part and in this way to reduce the dimension of the system. Furthermore, because 
the abelian differential Galois group implies in particular that this group is solvable, thus we can 
use directly all solvability results concerning some known equations such as e.g. hypergeometric 
equation, Whittaker equation. Lame equations. In addition, for a linear second order equation 
with rational coefficients there exists the closed algorithm, so-called Kovacic algorithm [4UJ, 
that decides whether equation is solvable in a class of Liouvillian function, yields explicit forms 
of solutions as well as determines the differential Galois group. This is achieved by providing 
necessary conditions for solvability of the appropriate Galois group. The equations in question 
have as their Galois group an algebraic subgroup of SL2(C), and since there are only three 
possibilities of those having a solvable identity component, the procedure is arranged in three 
cases. They consist of analysing the equation's singular points and finding an appropriate 
polynomial and an algebraic function of possible degrees 2, 4, 6, or 12; used to construct solutions. 

This means that Theorem [1] yields really an effective tool for proving the non-integrability 
and distinguishing the cases suspected about integrability in the case when Hamiltonian depends 
on some physical parameters, for examples of applications see references in [52]. 




It can happen that a considered system satisfies all conditions of the above theorem, but 
nevertheless it is not integrable. It is nothing strange as this theorem gives only necessary 
conditions for the integrability. This shows a need of stronger necessary conditions for the 
integrability. They were developed by C. Simo, J.J Morales and J.-P. Ramis [l9l[50l[52] and are 
based on higher order variational equations (HVE's). 

Theorem 3 Assume that a Hamiltonian system is meromorphically integrable in the Liouville 
sense in a neighbourhood of the analytic phase curve T, and the infinity is a regular singular 
point of the variational equations along T. Then the identity component of the differential Galois 
group of k-th variational equations along T is Abelian for all k > 1. 

For variational equations (VE's) of degree greater than one there is no more the splitting of 
variational equations into two parts: normal (NVE's) and tangential (TVE's) and there is no 
reduction of system's dimension and the analysis of the differential Galois group of the whole 
system of variational equations is more involved. Fortunately, in the case when variational 
equations are the product of Lame equations with Lame-Hermite solutions Morales- Ruiz proved 
in |49t [Sni E2] that the absence of logarithmic terms in solutions of higher order variational 
equations is a necessary condition of abelianity of the identity component of their differential 
Galois groups. 



The interested reader more detailed and complete presentation of Morales-Ramis theory can 
find in da [Ml ED ESI |66] and of differential Galois theory in H ES EZl SHl ES] . 
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4 Analysis of the minimally coupled field 



4.1 A = case 

The system now has the fohowing form 



a- 



-Pa, Pa = '2.ka - Arri^a^ct)^ + -^P% 
\ 

-P<t>, P(i> = -2m^a'^((>. 



(24) 



Using the aforementioned particular solution, for which the constant energy condition becomes 



we have as the variational equations 



/ a(^) \ / -1 
2k 



Pa 







V -2m'^q'^ / 



\ / a^'^ \ 


-2 



Pa 



(25) 



The normal part of the above system, after eliminating the momentum variation p^^"* , and writing 
for (^(1), IS 



qx + 2qx + 2m^q^x = 0, 



(26) 



which we further simplify like before by taking z = q as the new independent variable, and using 
the energy condition to get 



z{E - kz'^)x" + {2E - 3kz^)x' + m^z'^x = 0. 



2\ / 



(27) 



We check the physical hypersurface oi E = 0. This requires k ^ for otherwise the special 
solution would become an equilibrium point. Introducing a new pair of variables 



w{s) = w ( ^y^^ 



m 



,3/2 



x{z), 



we finally get 



l_K 4//^ - 
ds2 ~ V 4 s 4s2 



w, 



(28) 



(29) 



with fj, = ±1, and k, = 0. This is the Whittaker equation, and its solutions are Liouvillian if, and 
only if, («; + ;U — |, K — /X — ^) are integers, one of them being positive and the other negative 
[49j . As this is not the case here, this finishes the proof for k ^ 0. We recall, that because of 
the irregular singular point s = oo, this rules out only the rational first integrals. 

Non-integrability on one energy level means no global integrability, for the existence of 
another integral for all values of E would imply its existence on ii^ = 0. However, there might 
exist additional integrals for only some, special values of the energy. It is straightforward to 
check with the use of Kovacic's algorithm [3{)J, that this is not true here. For our equation, in 
cases 1 and 2 of the algorithm, there is no appropriate integer degree of a polynomial needed for 
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the construction of the solution, and case 3 cannot hold, because of the orders of the singular 
points of the equation. 

If /c = 0, a change of the dependent variable to w{z) = zx{z), reduces equation (f271) to 



Ew" + m^z^w = 0, 



(30) 



which is known not to possess Liouvillian solutions 

We notice that when A = = A; = 0, the system can be reduced to a two-dimensional one. 
In fact, the reduction is still possible when A 7^ 0, so we choose to present in the next section. 



4.2 A 7^ case 

We use the nonzero constant A to rescale the system as follows 

_ gi _ Pi 

""^/A' ^'^"VA' 



92, 



P2 

A'' 



(31) 



so that the equations become 



qi = -Pi, 



pi=2kqi-4qf{l+bq^) + \pl 

qf 



(32) 



Q2 = -^P2, 

qr 



P2 



-2bq2qt, 



where b = m^/A. The energy integral, for the previously defined particular solution, now reads 
£ = EA = ■^q'^ + kq^ — g^, where q has been rescaled according to (f3T]) . 
As before, we are interested in the variational equations, which read 



/ q? \ 

P? 
qi'' 
\ P^' J 



( -1 

2(k - 6g2) 











-2bq^ 



\ 


/ 



\ P? J 



(33) 



and writing x for q^\ and y for ^2^^ The normal part is 

1 

q^ 

ij = —2bq^x, 



(34) 



or alternatively 



X + 2-x + 2bq^x = 0. 



(35) 
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4.2.1 £ = 



We first pick the particular solution lying on the zero-energy level, as the global integra- 
bility implies the integrability for this particular value of the Hamiltonian. It is important to 
remember, however, that the converse is not true. 

The normal variational equation is cast into a rational form by changing the independent 
variable to z = q'^/k (for k ^ which implies k'^ = 1), and using the energy first integral. It 
then becomes a hypergeometric equation 

„ 5z — 4: , b 
^+2^(^^+41(7^^ = °' 

with the respective characteristic exponents 

z = 0, p=-l,0 

^ = 1' P = ^^\ (37) 
z = oo, p= ^(3-\/9^^),^(3 + \/9^^). 

By Kimura's theorem [3S|, the solutions of equation (j36|) are Liouvillian if, and only if 9 — 46 = 
(2p — 1)^, p e Z. As before, this means that for the global integrability this condition must be 
satisfied. 

For k = the solution of NVE is xi^2 = q~'^P°°^-^ , and the reduction to a two-dimensional 
system is possible, as mentioned before. 

4.2.2 

The special solution, is now directly connected to the Weierstrass p function, for if we 
introduce a new dependent variable v with 

'^' = r + ^' ^^^^ 
the energy integral implies that it satisfies the equation 

= 4i;^ - g2V - 93, (39) 

where 

92 = fik'- 93 = ^k{2e - 9£), (40) 

and the discriminant A = W2A£'^{k'^ — 4iS), which we take as non-zero to consider the generic 
case, 
reads 



case. Thus, taking w = q^\, and eliminating p^^ as before, the normal variational equation 



w = [Ap{'n;92,93) + B]w, (41) 

with A = 2 — b and B = — |A;(1 -|- b). This is the Lame differential equation, whose Liouvil- 
lian solutions are known to fall into three mutually exclusive cases, which are exactly those of 
Kovacic's algorithm: 
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1. The Lame-Hermite case, with A = n{n + 1) = 2 — 6, n G N. This imphes that 9 — 46 = 
(2n + 1)^. Case with n = 1 aheady known to be integrable because 6 = represents the 
massless field. 

2. The Brioschi-Halphen-Crawford case, where necessarily n is half an integer, i.e. n + ^ = 
I € N, and as before 9 - 46 = (2n + l)^ = {2lf. 



3. The Baldassarri case, with n + | S |Z U |Z U 
on B,g2, and 53. 



1 \ Z, and additional algebraic restrictions 



In the Lame-Hermite case we have infinite number of values of 6 = 2 — n(n + 1), n G N, for 
which the necessary conditions for the integrability given by the Morales-Ramis Theorem [T] 
are satisfied. In order to obtain stronger result we need to apply a more restrictive necessary 
conditions. Such conditions are given by higher order variational equations, see for detailed 
exposition. Here we explain this technique on the considered problem and we will follow |45j . 

At the beginning, it is convenient to change the variables in equations (j32p in the following 
way 

qi = wi, pi = -W2, 



W3 

92 = — , P2 



WIW4 — W2WS. 



Let 



w 



W{w), W = {wi,W2, 11)3,104^), 



(42) 



(43) 



be the system ()32p written in the new variables. The advantage of new coordinates is that now 
the variational equations split into a direct product of two Lame equations 



w 



w 



(1) 

1 

(1) 

2 

(1) 



V / 



/ 1 \ 

Aip{r])+Bi 

1 

V A2p{v)+B2 / 



(1) 
1 
(1) 



w 



w 



(1) 

3 

(1) 



(44) 



where p(r/) is the one given by equation ([HUj) . and 



^1 = 6, 
A2 



Bi = 2k, 



n(n+l), B2 = -kin"^ + n-3). 
o 

To derive the higher order VE's we substitute into equation ()43|) the infinite formal series 



(45) 



2,„(2) , ,3,„(3) 



(46) 



where 99 is the particular solution, and get 



1,2, 



(47) 
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where W{(p{r])) is the matrix of right hand sides in (I44p . and fj^w 



(1) 



,w 



(i-i) 



) are vectors 



obtained from the Taylor expansions of components of W{w). In particular we have 
/i = 0, 

b 



(48) 



and so on. For j = 1 equation (|44|) is recovered. Although tD^^^, . . . , w^^^'^'^ enter polynomially 
in the right hand sides of j-ih. variational equations (I47p , there exists an appropriate framework 
to define their differential Galois group. In [52] it was proved that if the system is integrable, 
then the identity component G° of the differential Galois group Gj of j-th variational equations 
is Abelian. Generally it is very difficult to determine Gj for j > 1. However, in a case when the 
first variational equations are a product of two Lame equations having infinite differential Galois 
group we have an effective method to decide whether is Abelian. Namely, if a logarithmic 
therm appears in local solution around r/ = of j-ih. variational equations, then is not 
Abelian, see [52l [53] for details. 

The calculations proceed as follows. The solution of ([T7|) is given by 



w 



iJ) 



(49) 



where fj = fj{w^^\ . . . ,w^^ ^^) and X is the fundamental matrix of the homogeneous system 
(i.e. the first order VE (dH)), so that 



X = W'{^{t))X, detX/0. 



We took 



X 



( Vx V2 
Vi V2 



\ 




Vs V4 
V V3 V4 J 



(50) 



(51) 



with 

Vl 
V3 



on + 9 



V2 



1 k_ 

brf' 15 
1 



+ 



k{v? + n - 3) 



(2n + l)r/" (2n + l)(6n - 3)7?"-2 



+ 



Next, we take as the solution of the first order VE 

'(1) = (0,0,t;4,t'4). 



(52) 



(53) 



Then we fix n = 2 and solve the second order VE and we obtain the integrand of (j49p for j = 3 
to be 

1 . . 54 44A; 54 128A; , , 

X-V3 = (0,0,— -^ + -- - ,—_-__ + ...), (54) 



625r/« 625??^ 



' 1257?3 875?? 
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which produces a logarithm in w^^\ If k = 0, one has to find solutions of fourth order VE to get 



1 r _ (f. 3618 1272g 3618 1536g 

h- (U, U, ^Qg375^io 218757?6 + ' ' ' ' 21875?75 218757? ^ ' 

which proves the non-integr ability, since we assumed £ ^ 0. 

This behaviour does not change as we increase n, although it was checked only for 10 con- 
secutive values. We thus conjecture that for 6 = 2 — n{n + 1) with integer n > 1 the system is 
not integrable. The procedure described is correct under assumption that the differential Galois 
group of the normal variational equations is not finite. We discuss this point in Appendix C, 
and justify that except countable many values of energy the group is not finite. 

In the Brioschi case, there is another additional condition for the integr ability: the so-called 
Brioschi determinant Qi is zero [39] . Unfortunately, there is no closed formula for Qi for general 
/, but analysing the first few values we notice a pattern: 

Qi = -^k, 

Q3 = -lk{35k'^ - 192S), 

\ (56) 

Q4 = (2835/c^ - 21600/c2f - 48384£:2), 

16 

A TO ^ 

Q5 = ^A;(231A:^ - 2240^^^ - 16384^2)^ 

8505 

Qe = — (15015/c^ - 176400fc''£: - 2802432A;^£:2 - 1126400£'^). 

64 

When k = 0, Qi is zero for odd /, and proportional to energy, which is not zero, for even I. When 
A; 7^ 0, so that k'^ = 1, each Qi is a polynomial in £, and that gives at most a finite number of 
energy values for which Qi = and the system is potentially integrable. We, again, conjecture 
that if the system is integrable (with this subsection's assumptions) and A; = 0, then necessarily 
n + ^ is odd, and that if /c^ = 1, then it is not integrable on a generic energy level. 

The Baldassarri case can also be studied in more detail by means of the modular function 

_ gl _ 4ie-3£f forP = l, 

^ gl-Tlgl 21{k^-A£)£^ \i for A = 0. ^ ^ 

A theorem by Dwork [49j states that the number of pairs (j, B) is at most finite in integrable 
cases. Since j depends on the energy for non-zero k, and B depends on m?, it means that a 
generic energy level is not integrable for a given value of m?. 

4.2.3 8 = k = {) 

As mentioned in Section 14.11 in this case we can transform the system to a two-dimensional 
one. In order to do that, time needs to be changed from the conformal to the cosmological one 
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dr/ — > dt = adr] in the original equations ()lip . We then take as the new momenta the Bubble's 
function a and the derivative of (j) 

h ■- -— - 

'~adt~ a?' (5g) 
^ '~ dt ~ a^' 

(This uj is not to be confused with the one introduced in Section 2.) Accordingly we have 

ah, 



da 
dt 



dt 
dh 

'dt 

du) 

~dt 



(59) 

—Irn^cp — Sojh. 



Thus, we are left with a dynamical system in the {h, cp, uj) space, as a decouples. Furthermore, 
the energy integral is now 

= ^a^(2A + 6j2 + 2m^(/>2-/i2), (60) 

so for a(t) which is not trivially zero, it gives a first integral on the reduced space. Choosing an 
appropriate variable a, suggested by the form of this integral 

^ Vh^ - 2A . , , 

(p = — sm(a), 

V2m (61) 

UJ = \Jh^ - 2Acos(q), 

we finally obtain 

— = \/2 m + 3/i sin(a) cos(a), 

dl ^^'^ 
^ = _3(/i2_2A)cos2(a). 



The problem of such reduction was also discussed in [28]. It is argued that there can be no 
chaos in this system, but its integrability - which would be one more first integral - remains 
unresolved. 



5 Analysis of the conformally coupled field 

5.1 Known integrable families 

There are four known cases when the system has an additional first integral, functionally 
independent of the Hamiltonian. They were found by applying the so-called ARS algorithm 
basing on the Painleve analysis [1]. Table [1] summarises those results. 
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case 


k 


A 




(1) 


0,±1 


A = A 


= -3A 


(2) 


0,±1 


A = A 


m'^ = —A 


(3) 





A = 16A 


= -6A 


(4) 





A = 8A 


= -3A 



Table 1: Known integrable cases for the conformally coupled field 



And the respective integrals of the systems are 



(1) 



(2)^ 



(3) 



(4) 



PlP2 



3k), 



\iA-Pi) + \{(i2-Qi)-'^{4 



H 

I = qiP2 + q2Pi, 



H = \{pl- pI) - ^{iH - iHil + 
I = {q\P2 + q2Pi)P2 + -Q-qiqUqi - 2^1), 



H 



1 771 

i:{pI - pI) - -p7(89^ - 6gf g| + 
2^2 - 



4qiq2PiP2 + qhl - {ql - ^qDpl + -^^qliql -'^qlf 



(63) 



In this work, we will show, that the above are the only integrable cases, when m 7^ 0. An 
important point to note is that there is a complete symmetry with respect to interchanging A 
and A. It is a consequence of the fact, that there exists a canonical transformation of the form 



Pi^ipi, qi-^-iqi, 
P2^P2, q2^q2, 



(64) 



that changes the Hamiltonian into 
H 



\ [pI+pD + \ Wl + ql) - m\lq^\ + \ {Aqf + Ag|) , 



(65) 



which is the same after swapping the indices. We shall use this form of H, where the kinetic part 
is in the natural form, to make the use of some already existing theorems more straightforward. 



5.2 Integrability of the reduced problem 

It is possible to give stringent conditions for integrability of the system, by considering a 
reduced Hamiltonian. Namely, we can separate potential V into homogeneous parts of degree 2 
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and 4: 

V = Vh2 + Vm: 

Vh2 = \k {qf + qf) , ^66) 

VhA = \ {-2m\jq^2 + Wi + Mi) ■ 

The following fact is crucial in our considerations: if a potential V is integrable then its highest 
order as well as the lowest order parts are also integrable. This fact needs some additional 
justification as its several known proofs are not correct. In fact, consider potential V = Vmm + 
• • • + ^rnax, where Vyam and Fmax are homogeneous parts of ^ = ^(g), g G C" of the lowest and 
the highest degree, respectively. Assume that it admits meromorphic commuting independent 
first integrals Fi, . . . , Fn. If = Ri/Si for certain holomorphic functions Ri and Si, then we 
set fi = Ti/si, where and Sj are the lowest order terms of expansions of Ri and Si into the 
power series. It is easy to show that fi are first integrals of Vram- However, we cannot claim 
that they are functionally independent. Fortunately we can use in the described situation the 
Ziglin Lemma [65j which guarantees that we can always choose first integrals Fi in such a way 
that their lading terms fi are functionally independent. A more complicated situation arises 
with the integrability of Vmax- Here we have to assume that V is integrable with rational first 
integrals in order to distinguish their highest order terms. Then we need also an appropriate 
version of the Ziglin Lemma. Proofs of these facts will be published elsewhere. 

In our case if V given by ([66]) is integrable then Vh2 and Vm must also be integrable. Vh2 is 
the potential of the two-dimensional harmonic oscillator, thus, it is trivially integrable. However, 
the homogeneous part Vhi gives strong integrability restrictions for the whole potential V. We 
will call Vft4 the reduced potential and denote it by V . 

Thus we effectively set A; = 0, and are now in position to exercise known theorems concerning 
homogeneous potentials depending on two variables. In particular the complete analysis for 
degree 4 has been completed in [i5| . 

In order to identify our potential with some of the list given in that paper, we have to check 
how many Darboux points there exist, and what are the values of parameters A, A and m that 
give potentials equivalent to particular families. 

We say that a non-zero point (^i, 92) = d\s a, Darboux point of the potential V{qi,q2) when 
it satisfies the equation 

V'{d) = 7d, (67) 
where 7 G C* = C \ {0}. Such a point corresponds to a particular solution of the form 

q{rj) = f{r])d, p(r?) = /(r/)d, (68) 

with /(??) satisfying a differential equation that for a potential of degree 4 takes the form 

Rv) = -ifivf- (69) 

As explained in Section 3, particular solutions allow for studying the variational equations 
along them, and yield necessary conditions for existence of additional first integrals. However, 
the major simplification discovered in |45j is that additionally there is only a finite number of 
parameters' sets (or non-equivalent potentials) corresponding to integrable cases. 
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Following the cited paper's exposition and notation, we take 24.2 and I^^^ to be the sets of 
integrable homogeneous potentials of degree 4 with 2 and 3 simple Darboux points respectively. 
We recall also four characteristic potentials thereof 

V3 = \aqf + hqfq2 + ^{qj + qlf , 

V4 = \aq\ + 

4 ^ (70) 

V5 = + 3qUI + pi 

V, = 2q\ + Iqlql + \qi 

where a and h denote (for the sake of this paragraph) arbitrary complex numbers. 
We find that our potential has: 

1. Four simple Darboux points, when A(m^ + A)(m^ ~\~ ^) 7^ 0, and AA 7^ m^. The only 
integrable cases are: 

(a) A = A = —^iTL^ {V is equivalent to V4), 

(b) A = -|m,^, A = - \m? {V is equivalent to V5), 

(c) A = — A = — \rn? {V is equivalent to Vq). 

2. Three simple Darboux points, when A = 0, and A(m^ + '^) 7^ 0. There are no integrable 
families here as Z4^3 = 0. 

3. Two simple Darboux points, when either A = ^ and A(m^ + A) 7^ 0, or A = A = 0. 
Again, no integrable families are present here because 24^2 = 0- 

4. A triple Darboux point, when A = —m?. Additionally there is a simple Darboux point 
when A 7^ 0. The potential is equivalent to V^, and is only integrable when A = —m?. 

There are two immediate implications that follow. Firstly, that the main system itself with 
= is only integrable in those four cases, and the respective first integrals are known, as given 
in the table. Secondly, as it was shown in [33] those cases are the only ones which could be 
integrable when k ^ 0. This happens because the integrability of the full potential implies the 
integrability of the homogeneous parts of the maximal and minimal degree (the latter is trivially 
solvable in our case). 

As the table shows, when the potential is equivalent to V3 (or, to be precise, its integrable 
subcase) or V4, the second first integral is known; but V5 and Vq only have known first integrals 
with zero curvature. And as was shown in [12], for k = 1, the values of A and A are those of V5 
or Vg forbid integrability. This is easily extended to the k = —1 case, since after the change of 
variables 

g.^e-/4g., p. ^e-W4p., j = i^2, (71) 

we obtain a system with the sign of k changed, but the ratios m?/A and m^/A the same. Thus, 
concerning the conjecture of the quoted paper, our results for A; 7^ enable us to state, that it 
is true, when the rational integrability is considered. 



19 



However, the above considerations assume that the energy value is generic, so that the 
particular solution is a non-degenerate elliptic function. As stressed before, this does not preclude 
the existence of an additional first integral on the physically crucial zero-energy level. 



5.3 Integrability on the zero-energy level 

We choose not to investigate the Darboux points, but the variational equations directly, as 
they are considerably simpler in this case. The Hamiltonian equations of (j20p are 



(72) 



qi=Pi, pi = -kqi +m qiq2 - Aq^, 
q2=P2, P2 = -kq2 + mV,q2-Xql 

and they admit three invariant planes as was shown in [13]. They are 

Ilfc = {{qi,q2,Pi,P2) e C^lgfc = Apfc = 0}, A; = 1,2, 

Hs = UQ'i, 92,^1,^2) e C \q2 = aqi Ap2 = -api}, a = ^ , 

-|- A 

Obviously two particular solutions are 

{qi =Pi=0, q2 = q2{v), P2 = 92(^)}> ^ = ^ (^pj + kql + , 
{q2 =P2 = 0, qi= qiii]), pi = q[{r])}, = ^ (^+pj + kqf + ^q{^ , 



(74) 



and in order to find the third particular solution we make a canonical change of variables 

{qi,q2,Pi,P2f = HQi,Q2, Pi,P2f, (75) 

where symplectic matrix B has the block structure 

" A^f-? -:), 0^f° »V (76, 



A^r \-a b )' \0 

and a and b are defined by 



+ A , m? + X 



2m2-hA + A' y2m^ + X + A' 



(77) 



Let us introduce five quantities 



ai=2m^ + A + A, 02 = 3AA + 2m^(A + A) + m"^, 03 = v^A + m2)(A + m2) 
04 = A^ A^ - AA - m*, as = AA - m'^. 

Then, in the new variables, Hamiltonian (I20p has the form 



(78) 



H = ^ [Pf + Pi + k{Ql + QD] + ^ [a^Q\ + 2a2QlQl + 4(A - \)a^QiQl + a^Q^] , (79) 
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and the Hamiltonian equations read 

Qi = Pu Pi = -kQi - — [asQ? + asQiQi + (A - X)a3Ql] , 

7 (80) 

Q2 = P2, P2 = -kQ2 [a2QlQ2 + 3(A - X)a3QiQl + ^4^1] ■ 

ai 

Thus, the third particular solution can be seen to be 

{Q2 = P2 = 0, Qi = Qi(r/), Pi = Q[{v)}, = ^ (p? + kQl + ^Q^^ . (81) 

Of course, this is only valid for ai ^ 0. We investigate what happens when A + A = —2m? at 
the end of this section. 

Normal variational equations (NVE's) along those three solutions (in the position variables) 

are 

^"{r])=[-k + m^q{r]f] ^v), 
r(r?)= [-fc + mV^)2] Civ), 



(82) 



e(^), 



where q{ri) is one of {qiirfj, q2{'n), Qi(^)}) depending on the respective particular solution. 

We will consider the A; = case first. Changing the independent variable to 2; = q{r))^, all 
the NVE's are reduced to the following 

2z^eiz) + 3zC'{z)-Xiaz)=0, (83) 

whose solution is 

C{z) = ^-(i±VT+8A7)/4^ (g4^ 
where we have introduced three important quantities 

, , a2 3 - 2(Ai + A2) + AiA2 

Ai = T-, A2 = A3 = — = — . (85j 

A A as 1 — A1A2 

Note, that if any of A, A or 05 is zero, the corresponding particular solution is constant and 
cannot be used to restrict the problem's integrability. Thus, we are left with the E = k = case 
as suspected to be integrable. 

When wc assume 7^ 0, or more precisely k'^ = 1, and introduce the same independent 
variable z as before, the NVE's read 

2z'^{Az + 2k)C"{z) + z{3Az + 4fc)^'(z) + {m'^z - k)C{z) = 0, 

2z\Xz + 2k)e'{z) + z{3Xz + Ak)C'{z) + {m^z - k)C{z) = 0, ^gg^ 
2^2 (^^z + 2k^ i"{z) + z (^^z + 4fe^ i'{z) - [^z + k^ C{z) = 0. 

First, let us observe that unlike in the previous case, when any of A, A or is zero, the system 
is not integrable. This happens, because then the NVE's becomes the Bessel equation 

sT(s) + sC'(s) + (s'-n2)e(s) = 0, (87) 
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with n = 1 and in a new variable s = m\J zjk (for the first two equations) or s = m\f—2zfk (for 
the third equation). The Bessel equation is known not to possess Liouvilhan solutions for n = \ 
[40j . Together with the results of the previous section this leads us to the following lemma. 



Lemma 1 System h2U\) considered on the zero or generic energy level with k'^ = 1 is not inte- 
grable when A or X is zero. Additionally for A + A 7^ —2m?, it is not integrable when AA = m^. 

Assuming that none of those constants is zero, we rescale the variable z in the three equations 
with 

2k 2k 2kai , , 

z^ — -z, z^ — -z, z^ z, (88) 

A A 05 

respectively, so that all three are transformed into a Riemann P equation of the form 

^ l-7-7'\.,, , 77' /3/3'-M'-77'l , , , , 

+ + + — + r^^2+^^^^ TY^ =0, 89 

\ z z — I J [ (z — 1)^ z[z — I) J 

with the following pairs of exponents {6,5'), (7,7')) at its singular points 



(90) 



Using Kimura's results on solvability of the Riemann P equation [38j we check when the difference 
of the exponents give us cases with the necessary conditions for integrability satisfied, and find 
that the parameters must belong to the following families 

^_k{ii + il^ /.ez, i = 1,2,3. (91) 

These polynomials in li are invariant with respect to the change — 1, so it is enough to 

consider non-negative values only. Furthermore, Ai and A2 cannot be equal to zero, as m? ^ 0, 
so li and I2 need to be strictly positive. 

This result can be refined still. First, let us notice, that Aj are not independent and the 
relation between them is 

112 

+ ^ T + ^ T = (92) 



Ai - 1 A2 - 1 A3 - 1 

provided qi 7^ and 05 7^ 0. In the above form we had to exclude the possibility of Aj = 1, so 
we consider it separately. 

Both of Ai and A2 cannot be equal to 1, as that would mean 05 = and we have shown 
that then the equations are non-integrable if additionally ai 7^ 0. The ai = case is described 
below. 

If only one of Aj, say Ai is 1, then necessarily A3 = 1, which follows from the definition ()85p . 
and the only possibly integrable cases are those with A2 satisfying (j92p with I2 ^ 2. The same 
holds when Ai and A2 are interchanged. Also, A3 = 1 requires that one of the remaining Aj is 1. 

When li and I2 are taken to be grater than 1, Ai and A2 are positive, so the relation (|92|) 
requires that 2/(A3 — 1) is negative. This only happens for /3 = and it follows that h = I2 = 2, 
which is exactly the first known integrable case. Since l/(Ai — 1) and 1/(A2 — 1) are positive 
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and tend to zero monotonically as /j > 2 tends to infinity, there are no other solutions, and no 
other integrable sets of parameters. 

Finahy, we turn to see what happens when ai = 0, i.e. A + A = —2m?. This is equivalent to 

- + - = 2, (93) 

Al A2 



provided A 7^ and A 7^ and the same two conditions of ()9ip hold because the first two 
variational equations can still be used. It is straightforward to check that the only integer 
solution of 

1 1 _ 

/i(/i + l) ^ hih + l) " ^ ' 
is li = I2 = 1 (so, incidentally, 05 = 0), which we recognise as the second case of our table. 

6 Conclusions 

The main results of our paper can be summarised as follows. 
For the minimally coupled scalar fields, given by the Hamiltonian 

H = ^ + ^P^) - kqj + Aqf + w?qlq\, (95) 

we have: 

Theorem 4 For A = 0, if the system is integrable then necessarily E = k = 0. 

Theorem 5 When A 7^ 0, if the system is integrable on a generic energy level then either 

1. 9 - 4m^/A = for some I el^, or 

2. k = and9- Am'^/A = (2n + 1)^ /or n + 1 e IZ U IZ U IZ \ Z. 

Conjecture 5.1 Suppose A 7^ 0, and let n be an integer satisfying 9 — 4m^/A = (2n + 1)^. If 

the system is integrable on a generic energy level E ^ 0, then either 

1. n = 1 or n = —2 (m = in both cases), or 

2. k = and 9 — 4m^/A = (2/)^ with I an odd integer, or 

3. A; = anrf n + i E IZ U U \ Z. 

Note that this is more restrictive than Theorem [5l as case 1 of this theorem admits more values 
of n than the conjecture's cases 1 and 2 put together. 

Theorem 6 For the zero energy level, and provided that A 7^ 0, if the system is integrable then 
either 

1. k = 0, or 
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2. 9- W/A = (2n + if, neZ. 



While for the conformally coupled scalar fields, given by the Hamiltonian 

H = \ {-pI +pI) + \ [K-ql + 4) + m\lql] + \ {Kq\ + Ag|) , (96) 

we have: 

Theorem 7 The system restricted to a generic energy level E ^ is integrahle if, and only if, 

1. k = 0, and its parameters belong to the four families listed in TableUi Otherwise there 
exists no additional, meromorphic integral. 

2. k"^ = 1, and its parameters belong to the first two families of TableUl Other than that, 
there exists no additional, rational first integral. 

The second part of the above theorem can be strengthened to meromorphic first integrals, 
although not for all values of the parameters, as described in |12j . 

Theorem 8 // the system restricted to the zero energy level is integrable, then either 

1. k = {}, or 

2. k"^ = 1 and its parameters belong to the first two families of Table [II or 

3. k^ = I and one of {Xi, X2} is equal to 1, and the other satisfies the condition ( fPij) with 
k > 2. 

Otherwise, the system is not meromorphically integrable. In particular this means, that for 
k'^ = 1 if at least one of A and A is zero, then the system is non-integrable. 

These are, however, only necessary and not sufficient conditions, so that the system might 
still prove not to be integrable at all. In particular, the numerical search for chaos suggests both 
the lack of global first integrals, and crucial differences in the behaviour of the system for real 
and imaginary values of the variables. This might be a clue, that the system might have first 
integrals which are not analytic, and thus not prolongable to the complex domain. A system 
with similar property was studies by the authors in [46j . 

The immediate physical consequences of the non-integrability is the non-existence of con- 
stants or motion (by definition) or, in other words, laws of conservation. This results not only 
in the complexity of evolution but also in a harder descriptive approach to a physical system 
which does not possess any global, well defined, preserved quantities like total charge or spin (in 
general - we have not considered such quantities in the present work). It is also obvious that 
direct integration, or obtaining the solutions in closed forms by means of elementary functions 
is out of question with non-integrable problems. 

Of course, depending on the properties of the first integrals, we might get quite different re- 
sults, and the requirement of meromorphicity or rationality is still very restrictive. As described 
in the introduction, this leaves open the question of existence of real-analytic first integrals. 
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Also we recall that physically the scale factor a cannot even assume negative values, and some 
authors argue that when cosmological (instead of conformal) time is used, the evolution is not, in 
essence, chaotic [TS]. Thus, we would like to stress that Liouvillian integrability is a mathemat- 
ical property of the system, and often the methods used to study it require the complexification 
of variables. This means that when restricted to the narrower, physical domain, the dynamics 
might be much simpler. And in particular we might be interested in a particular trajectory 
whose behaviour is far from generic. It is no surprise then, that the dynamics of our system 
when restricted to a > might appear regular. It should still be noted that the notion of chaos, 
although frequently associated with the integrability, has not yet been successfully conflated 
with it. And that a regular evolution is not necessarily integrable. 
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Appendix A. Massless minimal field 

When m = 0, the Hamilton-Jacobi equation for the main Hamiltonian ()10p will become 
separable, because it can be written as 

^" "Ha*) -2° (ftr) -'="* + A«' + ;^, (97) 

with the full generating function S = W — Etj. Assuming W = A{a) + F{(j)), equation (j97p can 
be solved with 



A{a) = y J2 (^Aa^ - ka"^ - E + do. 



(98) 



where J is a constant of integration. The first equation of motion can then be deduced from 



dW f do 

r] = I — 1^^^^^^=^^^^^^= = const, (99) 

'2 [Ka^ - ka? - E + j,) 



dE 
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which can be rewritten as 



=2^Aa^-fca2-ii; + ^). (100) 



Or, introducing a new variable v = a^, as 



1 \ 2 

^) =8{Av^ -kv"^ -Ev + J), (101) 



so that the general solution is 



1 A; 

a^ = ^^ = ^p(f?-%;52,53) + ^, (102) 

where 

g^ = ^k^ + 16 AE 
32 64 

and ryo is the constant of integration. Of course, for A = 0, equation (jlOip admits solutions in 
terms of circular functions. 

The equation for (j){ri) is the following 

dW f del) f da , 

— = / , + / , = const, (104) 



^ a^j2{Ka^-ka^-E + j,) 



which we simplify by using the just obtained solution for v{rf) to get 



const = ^2^^+ (105) 

As V is an elliptic function of order two, the second integral can be evaluated by means of the 
Weierstrass zeta and sigma functions to yield 



(106) 



where r]i^2, are the zeroes of v(r]), given by 

3p(r?i,2;52,93) = -2/c, (107) 

and the constant of integration can be determined from the boundary conditions on the field cf). 
The functions ( and a are defined as follows 

-C'{z) = p{z), liin(C(^)-i) =0, 

a'{z) ^, , a{z) ^ (108) 

-^T = C(^), lm^ = l. 
(y\z) 2->o z 

Again, for J = 0, the integrals in (jl04p reduce to simpler functions. 
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Appendix B. Massless conformal field 

For m = we can separately solve equations for each variable, so that we have 



El = -la'^ - Ika'^ + ^Aa^ 



2 2 



2 -, 1 (109) 

2'^ ' 2(fP ' 2""^ ' 4' 



with El + E2 = E being the total energy. The first of these is immediately solved, when we 
substitute vi = c? to get 

vl = 2kv\-^kv\-%EiVi, (110) 

whose solution is 

2 2k 

viiv) = jp{v-m; 92,93) + (111) 

with 771 the integration constant and 

52 = lfc2 + 4AEi, g^ = ^k^ + ^kAEi. (112) 

Of course, when A = the Weierstrass function p reduces to a trigonometric function. 
Similarly, for the other variable, we substitute V2 = (t>^ and obtain 

vl = -2\v^ - Akv^ + ^Eiv - 4w'^ , (113) 

whose solution is 

2 2k 
V2{r]) = -jp{v -V2; 94,95) - 1^, (114) 

where 

54 = + 4AE2, 95 = ^k^ + ^k\E2 + XW, (115) 

and r/2 is the integration constant. As before, for A = the solution degenerates to trigonometric 
functions. 



Appendix C. Lame equation in the Lame-Hermit e case 

Let us consider Lame equation 

^ = {n{n + l)p{t) + B)y, n G N, (116) 

where the Weierstrass function has two periods 2u}i and 2u2 which are independent over M. We 
denote its differential Galois group over C{p,p) by G. 
Function v = p{t) satisfies the differential equation 

= Av^ - g2V - 93 =: f{v) (117) 



27 



The algebraic form of Lame equation is obtain from ()116p by setting z = p{t) and it reads 




y - 



n{n + l)z + B 

m 



y = 0, n E N. 



(118) 



Let Gal be the differential Galois group over C{z) of this equation. 

As it was show in [20, Section 5], G = Gal H SL(2, G), and moreover it was also shown that 
G is finite iff Gal is finite. In [9, Corollary 3.4] it was proved that if Gal is finite, then Gal is 
a dihedral group Dm of order 2m, for a certain m > 3. In this case. 



This fact implies that if G is finite, then it is a cyclic group of order m for a certain m > 3, so 
there are two independent solutions yi and y2 of (11160 such that E C(p, p), for i = 1,2. 

Now, it is known that for given n E N, and m > 3 the number of linearly non-equivalent 
Lame equations (I118p with differential Galois Dm is finite, see [HI [21]. Nevertheless, for long 
time it was unclear if there exists a Lame equation (jllSp for which Gal is finite. This problem, 
among other things, was analysed by Baldassarri and Dwork, see [2l[3lll], but only a bound on 
m was found. Later, see [51 E], examples of Lame equations (IllSp with a finite differential Galois 
group were found. 

In practice, it is important to distinguishing parameters n, B, §2 and for which Gal is 
Dm with prescribed m > 3. However, as far as we know, such conditions is difficult to obtain. 
For n = 1 and m = 5 such conditions are given explicitly in [5] where it is conjectured that for 
arbitrary n and m they should have a polynomial form with respect to variables B, §2 and (73. 
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